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We consider the case when decoherence is due to the fluctuations of some classical variable or 
parameter of a system and not to its entanglement with the environment. Under few and quite 
general assumptions, we derive a model-independent formalism for this non-dissipative decoherence, 
and we apply it to explain the decoherence observed in some recent experiments in cavity QED and 
- on trapped ions. 
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. I. INTRODUCTION 

Decoherence is the rapid transformation of a pure linear superposition state into the corresponding statistical 
mixture 

*>■ |V^)=a|V^i>+/5|V^2)^p™. = |anVi)(^i| + l/3|V2)(V^2|; (1) 

O ' 

■ this process does not preserve the purity of the state, that is, Trp^j^ < 1, and therefore it has to be described in terms 
T-H ' of a non-unitary evolution. The most common approach is the so-called environment-induced decoherence which 
, is based on the consideration that it is extremely difficult to isolate perfectly a system from uncontrollable degrees 
of freedom (the "environment"). The non-unitary evolution of the system of interest is obtained by considering 
the interaction with these uncontrolled degrees of freedom and tracing over them. In this approach, decoherence is 
— ' caused by the entanglement of the two states of the superposition with two approximately orthogonal states of the 
' environment lEi) and |i?2) 

PU. 

4!. ■ ia\iJi)+f3\iJ2))(E)\Eo)=>\ij)TOT = a\^Pi)®\Ei)+(3\^2)(E)\E2) . (2) 

^ . 

^ , Tracing over the environment and using {Ei\E2) — 0, one gets 
^3 . 

T^T^env {WtOt{iP\} = Pmix , (3) 
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• ^ , where pmix is defined in Eq. (|^). The environment behaves as a measurement apparatus because the states \Ei) behave 
/\ ' as "pointer states" associated with in this way the environment acquires "information" on the system state and 
therefore decoherence is described as an irreversible flow of information from the system into the environment ^ . In 
this approach, the system energy is usually not conserved and the interaction with the environment also accounts for 
the irreversible thermalization of the system of interest. However this approach is inevitably model- dependent, because 
one has to assume a model Hamiltonian for the environment and the interaction between system and environment. This 
modelization, and therefore any quantitative prediction, becomes problematic whenever the environmental degrees of 
freedom responsible for decoherence are not easily recognizable. 

Decoherence is not always necessarily due to the entanglement with an environment, but it may be due, as well, 
to the fluctuations of some classical parameter or internal variable of the system. This kind of decoherence is present 
even in isolated systems, where environment-induced decoherence has to be neglected. In these cases the system 
energy is conserved, and one has a different form of decoherence, which we shall call "non-dissipative decoherence" . 
In such cases, every single experimental run is characterized by the usual unitary evolution generated by the system 
Hamiltonian. However, deflnite statistical prediction are obtained only repeating the experiment many times and this 
is when decoherence takes place, because each run corresponds to a different random value or stochastic realization of 
the fluctuating classical variable. The experimental results correspond therefore to an average over these fluctuations 
and they will describe in general an effective non-unitary evolution. 

In this paper we shall present a quite general theory of non-dissipative decoherence for isolated systems which can 
be applied for two different kinds of fluctuating variables or parameters: the case of a random evolution time and 
the case of a fluctuating Rabi frequency yielding a fluctuation of the Hamiltonian. In both cases one has random 
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phases e"'^"*/'' in the energy eigenstates basis that, once averaged over many experimental runs, lead to the decay 
of off-diagonal matrix elements of the density operator, while leaving the diagonal ones unchanged. 

The outline of the paper is as follows. In Section II we shall derive the theory under general assumptions, following 
closely the original derivation presented in 1^,1). In Section III we shall apply this theory in order to describe the 
decoherence effects observed in two cavity QED experiments performed in Paris, one describing Rabi oscillations 
associated with the resonant interaction between a Rydberg atom and a microwave cavity mode |5[| , and the second 
one a Ramsey interferometry experiment using a dispersive interaction between the cavity mode and the atom 
In Section IV we shall apply our approach to a Rabi oscillation experiment for trapped ions Q , and Section V is for 
concluding remarks. 



II. THE GENERAL FORMALISM 

The formalism describing non-dissipative decoherence of isolated systems has been derived in by considering 
the case of a system with random evolution time. The evolution time may be random because of the finite time needed 
to prepare the initial state of the system, because of the randomness of the detection time, as well as many other 
reasons. For example, in cavity QED experiments, the evolution time is the interaction time, which is determined by 
the time of flight of the atoms within the cavity and this time can be random due to atomic velocity dispersion. 

In these cases, the experimental observations are not described by the usual density matrix of the whole system 
p(t), but by its time averaged counterpart ||,^ 

p{t)= r dt'P{t,t')p{t') , (4) 



where p(t') = exp{—iLt'}p{Q) is the usual unitarily evolved density operator from the initial state and L... = 
[H, . . .]/'h. Therefore t' denotes the random evolution time, while i is a parameter describing the usual "clock" time. 
Using Eq. (^), one can write 

m = v{t)p{o) , (5) 

where 

V{t) ^ / dt'P{t,t')e~'^*' (6) 



is the evolution operator for the averaged state of the system. Following Ref. we determine the function P{t, t') 
by imposing the following plausible conditions: i) p{t) must be a density operator, i.e. it must be self-adjoint, positive- 
definite, and with unit-trace. This leads to the condition that P{t, t') must be non-negative and normalized, i.e a 
probability density in t' so that Eq. is a completely positive mapping, ii) V{t) satisfies the semigroup property 
V{ti + ts) = Viti)V{t2), with ii, i2 > 0. 

The semigroup condition is satisfied by an exponential dependence on t 

V{t) = {Fi}-*/^^ , (7) 

where T2 naturally appears as a scaling time. A solution satisfying all the conditions we have imposed can be found by 
separating Vi in its hermitian and antihermitian part Vi = A + iB and by considering the Gamma function integral 
identity [g| 

(l/i)-*/-2 ^ (A + iB)-*/"^ = -^-^ / dW^-^e-^^e-'^^ . (8) 

r(^ Wo 



Now the right hand side of Eq. (||) can be identified with the right hand side of Eq. (|^) if we impose the following 
conditions: A = t'/ri, where ri is another scaling time, generally different from T2; B = Lti in order to make the 
exponential terms identical, and A — 1 in order to get a normalized probability distribution P{t,t'). This choice 
yields the following expressions for the evolution operator for the averaged density matrix V{t) and for the probability 
density P(t,i',ri,r2) |,| 

Vit) = {l + zLn)-'/^' (9) 
P{t,t',n,r2) = ^ V , ^ • (10) 
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Notice that the ordinary quantum evolution is recovered when ri = r2 = r ^ 0; in this hmit P{t, t', ti, T2) S{t — t') 
so that pit) = pit) and V{t) = exp{— iLt} is the usual unitary evolution. Moreover, it can be seen that Eq. (|) 
implies that p{t) satisfies a finite difference equation |^ . The semigroup condition leads to the form of the probability 
distribution i', ti, T2) we use to perform the average on the fluctuating evolution times. However, notice that 
this probability distribution depends on both the two scaling times ti and T2 only apparently. In fact, if we change 
variable in the time integral, t" — {T2/Ti)t\ it is possible to rewrite the integral expression for V{t) in the following 
way 

/"CXD 

V{t) = (1 + iLti)-''^^ = I dt"P{t, t'\ r2)e-'^("i/"^)*" , (11) 

where 







p{t,t'\r2) ^ - — , , . (12) 

T2 r(t/T2) 

This probability density depends only on T2. However Eq. (^) contains an effective rescaled time evolution generator 
Leff = L(ti/t2). The physical meaning of the probability distribution of Eq. (p^, of the rescaled evolution operator. 



and of the two scaling times can be understood if we consider the following simple example. Let us consider a system 
with Hamiltonian H{t) — f{t)Ho, where 

00 

fit) = ^ ^(t - nT2)e{nT2 +n-t) (13) 
n=0 

{9{t) is the Heaviside step function), that is, a system with Hamiltonian Hq which is periodically applied for a time 
Ti, with time period T2 {t2 > ti) and which is "turned off" otherwise. The unitary evolution operator for this system 
is U{t) = e-^-^(*)^o, where Lq ^ [Hq, . . .] and 

p/ .N _ / ^ + "('^1 -'^2) nT2 <t < nT2 + Tl , . 

(n+l)ri nT2 + Ti <t < {n+l)T2 ' ^ ' 

which can be however well approximated by the "rescaled" evolution operator Ueff{t) — e^'^^°*^^'^^^'^^\ In fact, the 
maximum relative error in replacing F{t) with i(ri/r2) is {t2 —Ti)/t and becomes negligible at large times (see Fig. 1). 
This fact suggests to interpret the time average of Eq. ( pd] ) as an average over unitary evolutions generated by L, taking 
place randomly in time, with mean time width ti, and separated by a mean time interval T2. This interpretation is 
confirmed by the fact that when t = kT2, for integer k, the probability distribution P(t, i", T2) of Eq. ( p^ is a known 
statistical distribution giving the probability density that the waiting time for k independent events is t" when T2 is 
the mean time interval between two events. A particularly clear example of the random process in time implied by 
the above equations is provided by the micromaser |9| in which a microwave cavity is crossed by a beam of resonant 
atoms with mean injection rate R = I/T2, and a mean interaction time within the cavity corresponding to Ti. In the 
micromaser theory, the non unitary operator M describing the effective dynamics of the microwave mode during each 
atomic crossing replaces the evolution operator 6"*^"^^ Another example of interrupted evolution is provided by 
the experimental scheme proposed in [ pO[ for the quantum non-demolition (QND) measurement pd| | of the photon 
number in a high-Q cavity. In this proposal, the photon number is determined by measuring the phase shift induced 
on a train of Rydberg atoms sent through the microwave cavity with mean rate 1/t2, and interacting dispersively 
with the cavity mode. These two examples show that the two scaling times ti and T2 have not to be considered as 
new universal constants, but as two characteristic times of the system under study. 

However, in most cases, one does not have an interrupted evolution as in micromaser-like situations, but a standard, 
continuous evolution generated by an Hamiltonian H. In this case the "scaled" effective evolution operator has to 
coincide with the usual one, L, and this is possible only if ti — T2 — t. In this case r is simply the parameter 
characterizing the strength of the fluctuations of the random evolution time. This meaning of the parameter t in the 
case of equal scaling times is confirmed by the expressions of the mean and the variance of the probability distribution 
of Eq. (0) 

{t') = -t (15) 

T2 

a^{t') ^ - (t')2 ^ lit . (16) 

T2 
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When Ti = T2, the mean evolution time coincide with the "clock " time while the variance of the evolution time 
becomes <J^{t') = tr. In the rest of the paper we shall always consider the standard situation of an isolated system 
with Hamiltonian H , continuously evolving in time, and we shall always assume ti — T2 = t. 

When r = 0, V{t) — exp{— zLi} is the usual unitary evolution. For finite r, on the contrary, the evolution equation 
describes a decay of the off-diagonal matrix elements in the energy representation, whereas the diagonal matrix 
elements remain constant, i.e. the energy is still a constant of motion. In fact, in the energy eigenbasis, Eqs. (H) and 
(i) yield 

= + - H^^l^riynr. - e--.-e-".-p„,.(0) , (17) 

where LiJn,m = {En — Em)/'h and 

ln,m^^\0g{l+ul„jl) (18) 

i^n.m = — arctan(a;„,mTi) . (19) 

T2 

This means that, in general, the effect of the average over the fluctuating evolution time yields an exponential decay 
and a frequency shift LJn,m —> Vn^m of every term oscillating in time with frequency ujn.m- 

The phase diffusion aspects of the present approach can also be seen if the evolution equation of the averaged 
density matrix p{t) is considered. In fact, by differentiating with respect to time Eq. (|^) and using (^), one gets the 
following master equation for p{t) (we consider the case ti = T2 = r) 

^(i) = -i log (1 + zLr) p(0 ; (20) 
expanding the logarithm at second order in Lr, one obtains 



m = -l[H, m] - [H, [H, Pirn , (21) 



which is the well-known phase-destroying master equation [[12[ . Hence Eq. (|20| ) appears as a generalized phase- 
destroying master equation taking into account higher order terms in t. Notice, however, that the present approach 
is different from the usual master equation approach in the sense that it is model-independent and no perturbative 
and specific statistical assumptions are made. The solution of Eq. ( pl| ) gives an expression for Pn,m(t) similar to that 
of Eq. (0), but with O] 



T 

7n,™ = ^ (22) 

which are nonetheless obtained also as a first order expansion 'm ti — T2 — r of Eqs. (|l|) and (|l|). The opposite 
limit w„i,nT ^ 1 has been discussed in detail in Ref. ||]. 

Finally a comment concerning the form of the evolution operator for the averaged density matrix V{t) of Eq. (|^). At 
first sight it seems that V{t) is in general a multivalued function of the Liouvillian L, and that V{t) is uniquely defined 
only when t/T2 = k, k integer. However, this form for V{t) is a consequence of the time average over P{t,t' ,ti,T2) 
of Eq. (p^), which is a properly defined, non-negative probability distribution only if the algebraic definition of the 
power law function (i'/ri)'-*/'^^-'"^ is assumed. This means that in Eq. ^ one has to take the first determination of 
the power-law function and in this way V(t) is univocally defined. 



III. APPLICATION TO CAVITY QED EXPERIMENTS 



A first experimental situation in which the above formalism can be applied is the Rabi oscillation experiment of 
Ref. 1^ , in which the resonant interaction between a quantized mode in a high-Q microwave cavity (with annihilation 
operator a) and two circular Rydberg states (|e) and \g)) of a Rb atom has been studied. This interaction is well 
described by the usual Jaynes-Cummings model, which in the interaction picture reads 

H = nnR{\e){g\a+\g){eW) , (24) 
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where flu is the Rabi frequency. 

The Rabi oscillations describing the exchange of excitations between atom and cavity mode are studied by inject- 
ing the velocity-selected Rydberg atom, prepared in the excited state |e), in the high-Q cavity and measuring the 
population of the lower atomic level g, Pf.g{t), as a function of the interaction time which is varied by changing 
the Rydberg atom velocity. Different initial states of the cavity mode have been considered in ||^. We shall restrict 
only to the case of vacuum state induced Rabi oscillations, where the decoherence effect is particularly evident. The 
Hamiltonian evolution according to Eq. (pi) predicts in this case Rabi oscillations of the form 



Pegit) = T;{^- cos {2nRt)) . (25) 



Experimentally instead, damped oscillations are observed, which are well fitted by 

= cos (217^0) , (26) 

where the decay time fitting the experimental data is — 40/isec and the corresponding Rabi frequency is 
Q.r/2ti — 25 Khz (see Fig. 2). This decay of quantum coherence cannot be associated with photon leakage out of 
the cavity because the cavity relaxation time is larger (220 /isec) and also because in this case one would have an 
asymptotic limit P^gP{oo) = 1. Therefore decoherence in this case has certainly a non dissipative origin, and dark 
counts of the atomic detectors, dephasing collisions with background gas or stray magnetic fields within the cavity 
have been suggested as possible sources of the damped oscillations. [p|Jl4|]. 

The damped behavior of Eq. (|2^) can be easily obtained if one applies the formalism described above. In fact, from 
the linearity of Eq. (^, one has that the time averaging procedure is also valid for mean values and matrix elements 
of each subsystem. Therefore one has 

/•oo 

Peg{t) = / dt'P{t,t')P,g{t') . (27) 

Using Eqs. (||), (^, ( p^ and (p5|), Eq. ( p7|) can be rewritten in the same form of Eq. ( p6| ) 

Pegit) ^^{l-e--" cos (lyt)) , (28) 



where, using Eqs. ( |18[ ) and (|19|), 

j^^\og{l + mlr^) (29) 
1/ = iarctg(2f]fl.r) . (30) 

T 

If the characteristic time t is sufficiently small, i.e. ftjiT ^ 1, there is no phase shift, ly ~ 217^, and 

7 = 2r2flT (31) 

(see also Eqs. ( |2^ ) and (23)). The fact that in Ref. Q the Rabi oscillation frequency essentially coincides with the 
theoretically expected one, suggests that the time r characterizing the fluctuations of the interaction time is sufficiently 
small so that it is reasonable to use Eq. (|3^). Using the above values for 7 and flfi, one can derive an estimate for t, so 
to get T ~ 0.5 fxsec. This estimate is consistent with the assumption ^Irt <C 1 we have made, but, more importantly, 
it turns out to be comparable to the experimental value of the uncertainty in the interaction time. In fact, the 
fluctuations of the interaction time are mainly due to the experimental uncertainty of the atomic velocity v, that is 
dt/t ~ Sv/v = 1% (see Ref. [Q), and taking an average interaction time i ~ 50 /isec, one gets r ~ 5t = tSv/v = 0.5 
/Ltsec, which is just the estimate we have derived from the experimental values. This simple argument supports the 
interpretation that the decoherence observed in |^ is essentially due to the randomness of the interaction time. In 
fact, in our opinion, the other effects proposed as possible sources of decoherence, such as dark counts of the atomic 
detectors, dephasing collisions with background gas or stray magnetic fields within the cavity, would give an overall, 
time-independent, contrast reduction of the Rabi oscillations, different from the observed exponential decay. 

Results similar to that of Ref. Q have been very recently obtained by H. Walther group at the Max Planck Institut 
fiir Quantenoptik, in a Rabi oscillation experiment involving again a high-Q microwave cavity mode resonantly 
interacting with Rydberg atoms [E5|. In this case, three different initial Fock states |n) of the cavity mode, n = 0, 1, 2, 
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have been studied, and preliminary results show a good quantitative agreement of the experimental data with our 
theoretical approach based on the dispersion of the interaction times. 

Another cavity QED experiment in which the observed decay of quantum coherence can be, at least partially, 
explained with our formalism in terms of a random interaction time, is the Ramsey interferometry experiment of M. 
Brune et al. [||. In this experiment, a QND measurement of the mean photon number of a microwave cavity mode is 
obtained by measuring, in a Ramsey interferometry scheme, the dispersive light shifts produced on circular Rydberg 
states by a nonresonant microwave field. The experimental scheme in this case is similar to that of the Rabi oscillation 
experiment, with two main differences: i) two low-Q microwave cavities i?i and R2, which can be fed by a classical 
source S with frequency lor, are added just before and after the cavity of interest C; ii) the cavity mode is highly 



detuned from the atomic transition {5 
picture with respect to 



3> ^r), so to work in the dispersive regime. In the interaction 



TlLjR 



[\e){e\~\9){g\] 



(we use the classical field as reference for the atomic phases), the Hamiltonian has the following dispersive form ||T 



H 



[\e){e\-\g){g\] + hx{t){\g){e\ + \e){g\) 



{\g){gWa- \e){e\aa 



(32) 
(33) 



where A — ujeg — lor, the Rabi frequency within the classical cavities x(i) is nonzero only when the atom is in Ri 
and i?2) and QR{t) is nonzero only within C. In the experiment, single circular Rydberg atoms are sent through the 
apparatus initially prepared in the state |e), and let us assume that the microwave cavity mode in C is in a generic 
state J2n Cnl'iT')- The atom is subject to a 7r/2 pulse in so that 



Cn\n)\e) 



l5» 



V2 



(34) 



Then the atom crosses the cavity C with an interaction time tint and the dispersive interaction yields 

-iujntint 



E 



V2 



fti„t+i-f{n+l)ti„ 



(35) 



Finally the atom is subject to the second tt/2 pulse in the second Ramsey zone R2 and the joint state of the Rydberg 
atom and the cavity mode becomes 



= ^c„e-^""^+''^*'-|7i)|cos 



nl, 1, A 



1, A ■ 



\g) 



(36) 
(37) 



where T is the time of flight from i?i to i?2- The experimentally interesting quantity is the probability to find at the 
end the atom in the g state, Peg{n,T), whose theoretical expression according to Eq. (^7\j is 



Fe<,(n,T) =cos2 



|(A-e„) 



where the photon number-dependent frequency shift e„ is given by 

-(2n+ 1) 



S d(i?i,i?2) 



(38) 



(39) 



In Eqs. (|3^) and ( |39[) we have used the fact that Unt/T is equal to the ratio between the waist of the cavity mode w 
and the distance between the two Ramsey cavities c?(i?i, i?2)- The actual experiment of Ref. Q has been performed in 
the bad cavity limity T^ei < Unt in which the cavity C relaxation time T^ei is smaller than the atom-cavity interaction 
time. In this case, the cavity photon number randomly changes during tint and in the corresponding expression ( |39| ) 
for the frequency shift e„, the photon number n has to replaced by the mean value n. The Ramsey fringes are observed 
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by sweeping the frequency of the classical source lor around resonance, that is, studying Peg(n,T) as a function of 
the detuning A. The experimentally observed Ramsey fringes show a reduced contrast, which moreover decreases for 
increasing detunings A (see Fig. 2 of Ref. Therefore one can try to explain the reduced contrast, i.e., the loss 
of quantum coherence, in terms of a fluctuating evolution time, which in this case means a random time of flight T 
originated again by the dispersion of the atomic velocities. We average again the quantity Peg of Eq. ( ^8| ) over the 
probability distribution P{t, t') derived in Section II, replacing t' with a random time of flight T', and we obtain 

Pe,(A) = i {1 + F(A - e„) cos [(A' - e^) T]} (40) 
where, using Eq. (^7|), the fringe visibility function i^(A — e„) is given by 

F(A - e„) = [1 + (A - e J '^'^^ , (41) 

and (A' — e'^) is the frequency shift 

(A'-e;) = iarctg[(A-e„)T] . (42) 

T 

The parameter t characterizing the strength of the fluctuations of the time of flight T can be estimated with arguments 
similar to those considered for the Rabi oscillation experiment. Since 5T/T ~ 5v/v — 1.5% and T ~ 300 /xsec (see 
Ref. Q), one has r ~ 5T ~ 4.5 /isec. For the interesting range of detunings A, one has (A — e„) r ^ 1, so that one 
can neglect again the frequency shift (^) and approximate the fringe visibility function (^Tj) with a gaussian function, 
that is, 

P,g(A)-i|l+e-^^^^^cos[(A-e„)T]| . (43) 

This gaussian modulation of the Ramsey fringes with a width cta — (Tr)^^/^ ~ 27 Khz is consistent with the typical 
experimental Ramsey fringe signal (see Fig. 2 of Ref. Q), but it is not able to completely account for the observed 
modulation and contrast reduction of the fringes. This means that, contrary to the case of the Rabi oscillation 
experiment, in this case the role of other experimental imperfections such as random phases due to stray fields, 
imperfect 7r/2 pulses in Ri and i?2 and detection errors, is as relevant as that of the dispersion of atomic velocities 
and these other effects have to be taken into account to get an exhaustive explanation of the observed decoherence. 



IV. RABI OSCILLATION EXPERIMENTS IN TRAPPED IONS 

Another interesting Rabi oscillation experiment, performed on a different system, that is, a trapped ion |Q, has 
recently observed a decoherence effect which cannot be attributed to dissipation. In the trapped ion experiment of 
Ref. 0, the interaction between two internal states (| t) and | J,)) of a Be ion and the ccntcr-of-mass vibrations in 
the z direction, induced by two driving Raman lasers is studied. In the interaction picture with respect to the free 
vibrational and internal Hamiltonian, this interaction is described by the following Hamiltonian [ p^ 

H ^ nn\ T) (i I cxp {i [r/ (ae-^^^* + flte*--*) - 6t + dp]] + H.C. , (44) 

where a denotes the annihiliation operator for the vibrations along the z direction, loz is the corresponding frequency 
and 6 is the detuning between the internal transition and the frequency difference between the two Raman lasers. The 
Rabi frequency Vl is proportional to the two Raman laser intensities, and rj is the Lamb-Dicke parameter 0,0. When 
the two Raman lasers are tuned to the first blue sideband, i.e. 5 = lOz, Hamiltonian ( |4^ ) predicts Rabi oscillations 
between | J,, n) and | ],n+ \) {\n) is a vibrational Fock state) with a frequency 

a. = ^-r=T^L\{Tf) , (45) 

where L\ is the generalized Laguerre polynomial. These Rabi oscillations have been experimentally verified by 
preparing the initial state | J,, n), (with n ranging from to 16) and measuring the probability Pj,(i) as a function of 
the interaction time which is varied by changing the duration of the Raman laser pulses. Again, as in the cavity 
QED experiment of ||] , the experimental Rabi oscillations are damped and well fitted by ||7|,p^ 
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Pli^: = ^(1 + e"''"* COS {2nnt)) , (46) 



where the measured oscillation frequencies r2„ are in very good agreement with the theoretical prediction (^5|) cor- 
responding to the measured Lamb-Dicke parameter r/ = 0.202 [Q. As concerns the decay rates 7„, the experimental 
values are fitted in by 

7« =7o(ri + l)°-^ (47) 

where 70 = 11.9 Khz. This power-law scaling has attracted the interest of a number of authors and it has been 
investigated in Refs. [TtIJi^ , even if a clear explanation of this behavior of the decay rates is still lacking. On the 
contrary, the scaling law (p7|) can be simply accounted for in the previous formalism if we consider the small r limit of 



Eq. (31), which is again suggested by the fact that the experimental and theoretical predictions for the frequencies r2„ 
agree. In fact, the n-dependence of the theoretical prediction of Eq. ( |4^ ) for 77 = 0.202 is well approximated, within 
10 %, by the power law dependence (see Fig. 3) 

r!„~r!o(« + l)°-'' , (48) 
so that, using Eq. ([3l|) , one has immediately the power law dependence (n-|-l)"-^ of Eq. (B^ . The value of the parameter 



T can be obtained by matching the values corresponding to n = 0, and using Eq. (|3l|), that is t = 7o/2rig ~ 1.5 • 10"^ 
sec, where we have used the experimental value r2o/27r = 94 Khz. 

However, this value of the parameter r cannot be explained in terms of some interaction time uncertainty, such as 
the time jitter of the Raman laser pulses, which is experimentally found to be much smaller |p^ . In this case, instead, 
the observed decoherence can be attributed, as already suggested in ]l6|-p^, to the fluctuation of the Raman laser 
intensities, yielding a fluctuating Rabi frequency parameter VL{t) of the Hamiltonian (|4^). In this case, the evolution 
is driven by a fluctuating Hamiltonian H[t) = fiVL{t)H, where H = H/il in Eq. (Q), so that 

Pit) - exp{-zZ /* d^m}p{0) = e-'^-^(*V(0) (49) 
Jo 

where L = [H , . . and we have defined the positive dimensionless random variable A{t) = J^d^n{^), which 
is proportional to the pulse area. It is now easy to understand that the physical situation is analogous to that 
characterized by a random interaction time considered in the preceding sections, with L replaced by L and t' by A(t). 
It is therefore straightforward to adapt the formalism developed in Section II to this case, in which the fluctuating 
quantity is the pulse area A, yielding again random phases in the energy basis representation. In analogy with Eq. (|^), 
one considers an averaged density matrix 

p{t)^ dAP{t,A)e-'^^'^p{<d) . (50) 







Imposing again that p{t) must be a density operator and the semigroup property, one finds results analogous to 
Eqs. d) and (0) 

V{t) = (l + imr) (51) 



Here, the parameters f2 and r are introduced as scaling parameters, but they have a clear meaning, as it can be easily 
seen by considering the mean and the variance of the probability distribution of Eq. (BT 



{A) ^ nt (53) 

a^{A) = (A^) - = QHt (54) 

implying that fl has now to be meant as a mean Rabi frequency, and that t quantifies the strength of A fluctuations. 
It is interesting to note that these first two moments of P{t, A) determine the properties of the fluctuating Rabi 
frequency i^(t), which can be written as 

n(t) = n + ^{t) (55) 

m = o mm) = ^^rs{t-t'), (56) 
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that is, the Rabi frequency ^l{t) is a white, non-gaussian (due to the non-gaussian form of P{t, A)) stochastic process. 
In fact, the semigroup assumption we have made imphes a Markovian treatment in which the spectrum of the laser 
intensity fluctuations is flat in the relevant frequency range. This in particular implies that we are neglecting the 
dynamics at small times, of the order of the correlation time of the laser intensity fluctuations. 

The estimated value of t gives a reasonable estimate of the pulse area fluctuations, since it corresponds to a fractional 
error of the pulse area -y/ a'^{A)/(A) = \/t jt of 10% for a pulse duration of t = 1 /isec, and which is decreasing for 
increasing pulse durations. 

The present analysis shows many similarities with that of Ref. jlj] which also tries to explain the decay of the Rabi 
oscillations in the ion trap experiments of Q in terms of laser intensity fluctuations. The authors of Ref. in fact 
use a phase destroying master equation coinciding with the second-order expansion ( pT| ) of our generalized master 
equation of Eq. ( pO[ ) (see Eq. (16) of Ref. ||l^ with the identifications G H/Ti and F ^ r) and moreover derive the 
same numerical estimate for the pulse area fluctuation strength F <-> r. Despite this similarities, they do not recover 
the scaling (47|) of the decay rates 7„ only because they do not use the general expression of the Rabi frequency 
(^5|), (and which is well approximated by the power law (^8|)) but its Lamb-Dicke limit ri„ = rio(". + 1)°'^, which is 
valid only when 7] 1. There is however another, more fundamental, difference between our approach and that of 
Ref. They assume from the beginning that the laser intensity fluctuations have a white and gaussian character, 
while we make no a priori assumption on the statistical properties of the pulse area A. We derive these properties, i.e. 



the probability distribution (52), only from the semigroup condition, and it is interesting to note that this condition 
yields a gaussian probability distribution for the pulse area only as a limiting case. In fact, from Eq. ( p2| ) one can 
see that P{t, A) tends to become a gaussian with the same mean value Vlt and the same width fJ^rt only in the large 
time limit i/r 3> 1 

P(i, ^=L= exp j (57) 

The non-gaussian character of P{t, A) can be traced back to the fact that P{t, A) must be definite and normalized 
in the interval < ^ < oo and not in — oo < ^4 < -|-oo. Notice that at i = r, Eq. assumes the exponential 
form P{t,A) — e~^^^'^ /fir. Only at large times t the random variable A becomes the sum of many independent 
contributions and assumes the gaussian form. 

Due to the non-gaussian nature of the random variable A, we find that the more generally valid phase-destroying 
master equation is given by Eq. (^0|) (with L replaced by QL) . The predictions of Eq. (^0|) significantly depart from 
its second order expansion in Lt, Eq. (|2l|), corresponding to the gaussian limit, as soon as r becomes comparable 
with the typical timescale of the system under study, which, in the present case, is the inverse of the Rabi frequency. 

The present analysis of the Rabi oscillation experiment of Ref. M can be repeated for the very recent experiment 
with trapped ions performed in Innsbruck pO| , in which Rabi oscillations involving the vibrational levels and an 
optical quadrupole transition of a single '^^Ca/ ion have been observed. Damped oscillations corresponding to initial 
vibrational numbers n = and n = 1 are reported. From the data with n — 0, Qq/2t: — 21 Khz and 70 = 1 Khz, we 
get r ~ 70 /2ilQ ~ 3 • 10^^ sec and this estimate is consistent with attributing again the decoherence to the fluctuations 
of the Rabi frequency caused by laser intensity fluctuations. Moreover in this case, the experiment is performed in 



the Lamb-Dicke limit 77 <C 1, and therefore, using again Eq. (31), we expect, in this case, a linear scaling with the 
vibrational number, 7„ — 20^r ~ 7o(n -I- 1). 



V. CONCLUDING REMARKS 



Decoherence is not always necessarily due to the entanglement with an environment, but it may be due, as well, to 
the fluctuations of some classical parameter or internal variable of a system. This is a different form of decoherence, 
which is present even in isolated systems, and that we have called non-dissipative decoherence. In this paper we have 
presented a model-independent theory for non-dissipative decoherence, which can be applied in the case of a random 
evolution time or in the case of a fluctuating Hamiltonian. This approach proves to be a flexible tool, able to give a 
quantitative understanding of the decoherence caused by the fluctuations of classical quantities. In fact, in this paper 
we have given a simple and unified description of the decoherence phenomenon observed in recent Rabi oscillation 
experiments performed in a cavity QED configuration |^ and on a trapped ion [Q. In particular, this approach has 
allowed us to explain for the first time in simple terms, the power-law scaling of the coherence decay rates of Eq. (|47|), 
observed in the trapped ion experiment. 

The relevant aspect of the approach applied here, and introduced in Ref. is its model-independence. The 
formalism is in fact derived starting from few, very general assumptions: i) the average density matrix p{t) has all 
the usual properties of a density matrix; ii) the semigroup property for the time evolution generator V{t) for p{t). 



9 



With this respect, this approach seems to provide a very general description of non-dissipative decoherence, in which 
the random properties of the fluctuating classical variables are characterized by the two, system-dependent, time 
parameters ri and T2. As we have seen in section II, in the cases where one has a standard, continuous evolution, 
the two times coincide ti — T2 = r. Under ideal conditions of no fluctuating classical variable or parameter, one 
would have t — 0, and the usual unitary evolution of an isolated system in quantum mechanics would be recovered. 
However, the generality of the approach suggests in some way the possibility that the parameter r, even though 
system-dependent, might have a lower nonzero limit, which would be reached just in the case of no fluctuations of 
experimental origin. This would mean a completely new description of time in quantum mechanics. In fact, the 
evolution time of a system t' (and not the "clock" time t) would become an intrinsically random variable with a 
well defined probability distribution, without the difficulty of introducing an evolution time operator. In Ref. it is 
suggested a relation of the nonzero limit for r with the "energy-time" fi/2/S.E, where /S.E is the uncertainty in energy. 
This would give a precise meaning to the time-energy uncertainty relation because now r rules the width of the time 
distribution function. However, this "intrinsic assumption" is not necessarily implied by the formalism developed in 
and apphed, with a more pragmatic attitude, in the present paper. 
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FIG. 1. The function F(t) defined in Eq. (0) (full line) is plotted as a function of time (expressed in arbitrary units) and 
compared with its "linear approximation", the rescaled time tri/r2 (dashed line). The relative error between them is given by 
(ti ~ Ti)/t and is negligible at large times t. 
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FIG. 2. The Rabi oscillations of the transition probability P,^g{t) as a function of time, according to the fitting function of 
Eq. M. 
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FIG. 3. The ratio between Rabi frequencies f2„/f2o experimentally measured in Ref. |7|, and well fitted by the theoretical 
prediction of Eq. (^), is plotted as a function of the initial vibrational number n (squares), and compared with the power- law 
approximation of Eq. (^), (n + 1)°'^^ (triangles). 
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